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The interplay of improper structural ferroelectricity (trimerization) and geometrically frustrated 
spins in the hexagonal RMnOs rare-earth manganites lead to a network of topologically protected 
defects that create strongly coupled structural and magnetic vortices. This unusual vortex structure 
is responsible for a significant domain wall magnetoelectric effect and magnetization, M, neither 
of which, however, occurs in the bulk. Here we demonstrate from a combination of first-principles 
calculations, group-theoretic techniques, and spin models that this structural trimerization not only 
exists in a family of hexagonal RFeOs ferrites, but also induces a bulk M, leading to a bulk coupling 
of ferroelectric and ferromagnetic domains such that if the polarization, P, reverses so does M. 
We show how this trimer distortion mediates an trilinear-coupling of P, M, and antiferromagnetic 
order, which is a universal feature of the hexagonal ferrites. 



Two themes at the forefront of materials physics are 
the cross-coupling of distinct types of ferroic order [1- 
|4 and topological defects in systems with spontaneous 
broken symmetry [5 -8 . Common to both are a plethora 
of novel phenomena to understand, and new properties 
and functionalities to exploit for novel applications. Mul- 
tiferroics |^ are an ideal platform to realize both 
themes in a single material. In this regard, an exciting 
development is the discovery of a topologically protected 
vortex-domain structure in one of the most extensively 
studied class of multiferroics, the hexagonal (hexa) rare- 
earth manganites. Here, antiphase structural ('trimer') 
domains are clamped to ferroelectric domain walls (and 
vice versa) pTHT4] forming a 'clover-leaf domain pattern. 
Fig. [T^. These trimer domains have a particular phase 
relationship that lead to an appearance of structural vor- 
tices, which in turn induce magnetic vortices [TSl [16], 
strongly coupling antiferromagnetism and the polariza- 
tion at the domain wall. This leads to a significant do- 
main wall magnetoelectric effect and magnetization nei- 
ther of which, however, occurs in the bulk. 

Our focus here is on elucidating a remarkable family of 
materials, the hexa RFeOa ferrites. We show they display 
all the ingredients necessary to realize the vortex-domain 
physics plus an intrinsic, bulk magnetization, M, induced 
by the trimer distortion, 0%^^ leading to entirely new 
magnetoelectric phenomena absent in manganites. 

The crystal structure of the hexa manganites is formed 
by a two-dimensional network of corner shared trigo- 
nal MnOs bipyramids, separated by layers of R ions. 
Fig. [TJd. Geometric frustration of the strongly antiferro- 
magnetic, AFM, nearest neighbor Mn spins leads to 120° 
non-collinear AFM ordering (Tat ^ 80-120 K). The key 
to realizing the unusual topological defects and strong 
coupling at the domain walls is the improper nature of 
ferroelectricity in this class of materials. Here the polar- 
ization, P, which is stable in the paraelectric PGa/mmc 
structure, is induced by a zone-tripling structural distor- 
tion, [UHlH]. The latter, referred to as the trimer 
distortion, is associated with a 2-up/l-down buckling of 



the R-planes and tilting of the MnOs bypyramids. It is 
nonlinearly coupled to the polarization, 

^trimer -PQk3COs(3^) (1) 

the form of which implies that a nonzero trimer distor- 
tion induces a nonzero P. There are three distinct ^ do- 
mains (a, /3, and 7) corresponding to one of three possible 
2c axes, i.e., one of the 3 permutations of 2-up/l-down. 
Also there are two tilting directions, either towards (+) 
or away from (— ) the 2c axis, i.e., l-up/2-down or 2-up/l- 
down, respectively. This results in six PGscm structural 
domains. A consequence of the improper origin of ferro- 
electricity is that the sign of P depends on the direction 
of Qxg- This simple fact leads to the nontrivial domain 
structure of the hexa manganites, Fig[l^ pTHT3l [T6]. 

Recently thin films of RFeOa have been epitaxially sta- 
bilized in the hexa PGscm structure [I9l[20]. These hexa 
ferrites exhibit ferroelectricity above room temperature, 
but with conflicting results as to its origin |2TJ [22] . As 
to the magnetic structure, experiments on LuFeOs sug- 
gest a considerably high AFM ordering temperature, Tn 
= 440 K [2F. Additionally there is evidence of a second 
magnetic transition around ~ lOOK, at which M becomes 
nonzero [2TH24j , however, the significance of this or even 
if it is an intrinsic or bulk effect is not previously known. 

Our main results are as follows: 1) ferroelectricity in 
the hexa ferrites is of the improper structural type where 
the trimer distortion induces P, and therefore a similar 
topological domain structure exists as in the manganites, 
however, the difference in electronic structure leads to 2) 
Fe spins ordering above room temperature and to 3) an 
120° AFM structure, L, that allows weak-ferromagnetism 
(M||z) induced by Q^g- The combined significance of 
which 4) suggests a bulk coupling of trimer/ferroelectric, 
FE, and ferromagnetic, FM, domains such that if P re- 
verses 180° so does M. A universal feature of all hexa 
ferrites is that the trimer distortion induces both P and 
M, and mediates an intrinsic bulk coupling of P, M, 
and L. These results are summarized in Figs, [ij^-e. In 
the remainder of this paper we explain these results and 
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FIG. 1. Improper structural ferroelectricity inducing 
weak-ferromagnetism, M^, in hexagonal ferrites from 
first principles. Results are shown for LuFeOa but (a-d) 
hold for RFeOs and RMnOs, while (e) occurs only for ferrites. 
a, Six structural domains of the primary "trimer" Q^g distor- 
tion and the secondary polarization P where P ^ Qf^^cosS^. 
Here <^=^ ^ = and a counterclockwise rotation corre- 
sponds to domains differing by ^ = +7r/3. Arrows indicate 
the direction of the trimer distortions, only distortions around 
the 2c axis are shown for clarity, b, Crystal structure of the 
(left) QJ^g = —Qks (center) paraelectric Qks = 0, and 

(right) — '^Qks domains respectively, c, Energy 

d, P, and e, as a function of Qxg, allowing only ^ = or 
TT trimer domains. Notice that P =0 when Q^g — 0? suggest- 
ing that a proper ferroelectric mechanism in not likely (see 
Supplemental). The insets of d, and e, show Energy vs. P 
and M;z, respectively. 



discuss the novel implications for the topological domain 
structure. 

Route to elucidate cross coupling of bulk or- 
der. Our approach to elucidate the bulk cross-coupling 
between ferroelectricity and magnetism is to introduce a 
PGa/mmc paraelectric, PE, reference structure that has 
the same AFM configuration as the PGscm FE ground 
state. The minimal model required to describe the spin 
structure [25] is the following effective spin Hamiltonian 

H = J2 JiJ^i ■ Si +1] • S,- X S,- +^ S, • *i • Si (2) 

ij ij i 

where the J^j's are the symmetric exchange interactions, 
Dij^s are the Dzyaloshinskii-Moriya (DM) antisymmet- 
ric exchange vectors, and is the single-ion anisotropy 
(SI A) tensor. Using this model as a guide, we then cal- 
culate from first principles how the material specific pa- 
rameters of this model change on going from the PE ref- 
erence structure to the FE ground state as a function of 
the trimer distortion. This also will allow us to explain 
why manganites and ferrites have qualitatively different 
magnetic order. Additionally, we will introduce a sim- 
ple phenomenological model that captures the essential 
physics, explaining our first-principles results and pro- 
viding an intuitive rational for the spin model. 

Method. The first principles calculations were per- 
formed using the DFT+U method [26] with the PBE 
form of exchange correlation functional [27]. We con- 
sidered Lu 4/ states in the core and for TM 3d states 
we chose U = 4.5 eV and Jh = 0.95 eV. Structural 
relaxations, frozen phonon and electric polarization cal- 
culations were performed without the spin-orbit coupling 
(SOC) using the projected augmented plane-wave basis 
based method as implemented in the VASP [28l [29] . We 
used a 4x4x2 k-point mesh and a kinetic energy cut-off 
of 500 eV. The Hellman-Feynman forces were converged 
to 0.001 eV/A. The electronic and magnetic properties 
were studied in the presence of SOC. 

We additionally cross-validated the electronic and 
magnetic properties using the Full-potential Linear Aug- 
mented Plane Wave (FLAPW) method as implemented 
in WIEN2K code [30]. 

The magnetic ground state. In the PE structure 
the transition metal, TM, spins form ideal triangular 
planes, which become slightly distorted in the FE phase. 
If the in-plane nearest neighbor (nn) TM-TM exchange 
interactions are AFM in nature, geometric frustration 
would cause the spins to order in a 120° non-collinear 
pattern, as in the manganites. Considering the com- 
patibility of this spin pattern with that of the crystal 
symmetry, four principle 120° magnetic spin configura- 
tions, denoted Ai, A2, Bi and B2, are allowed. Fig. [2^. 
Among these four principle phases only A2 (magnetic 
space group PGac'm') allows a net M along z. These 
spin states are defined by two free parameters, xi 
X2, as described in Fig. [2^. 
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FIG. 2. The noncollinear magnetic and electronic 
structure from first principles, a, The 120° antiferro- 
magnetic (AFM) ordering patterns - arrows indicate the spin 
direction within the xy plane. For each bipyramid layer, a = 1 
(red arrows, darker triangles) or 2 (blue arrows, lighter trian- 
gles), we define a local AFM order parameter, = L^c^e*^". 
The values of the angles (xi,X2), shown for the ^ = 
trimer domain, are defined by specific TM ions from adja- 
cent planes connected by the 2c axis that remains in the dis- 
torted phase. Among the four principle magnetic structures 
only A2 and Bi allow the xy spins to uniformly cant out of 
plane resulting in a net layer magnetization, MJ, = M^z^ and 
to weak-ferromagnetism (wFM) or weak- ant iferromagnet ism 
(wAFM), respectively. All intermediate spin patterns can 
be described as a linear combination of the four principle 
magnetic structures. First-principles calculations of b, Non- 
collinear magnetic energies and c, density of states (coUinear). 
Note that the difference in orbital occupancy is the fundamen- 
tal reason why ferrites, top, and manganites, bottom, have 
different magnetic ground states and subsequently the key to 
the predicted bulk coupling of polarization and magnetism in 
ferrites. Insets: crystal field splitting and occupancy of the 
TM majority channel. 



In addition to the principle magnetic configurations 
we also considered the four known intermediate mag- 
netic structures, Fig. |2^. The results of our total energy 
calculations for LuFeOs, LFO, and LuMnOa, LMO, are 
presented in Fig. (for clarity we limit our discussion 
to these two compounds). In agreement with non-linear 
optical measurements [31], we find that LMO stabilizes 
in the Bi state. Interestingly, this state displays weak- 
antiferromagentism, wAFM, i.e., there is a small canting 
of the spins out of the x7/-plane, which leads to a net mo- 
ment in each layer but consecutive layers cant in opposite 
directions leading to a net = 0. Note, however, that 
the B2 state, where the net magnetic moment is equal to 
zero in each layer by symmetry, is also close in energy. 

In contrast, LFO stabilizes in the A2 state where a 
canting of the spins out of the x7/-plane is also allowed 
but here consecutive layers add, leading to a net = 
0.02/iB/Fe along the z-axis, i.e., LFO displays weak- 
ferromagnetism, wFM [32^2- As shown in Fig.[2]3, how- 
ever, the Ai (M^ = by symmetry) and A' (finite al- 
lowed) states, lie energetically very close to ground state. 
It was recently inferred from neutron diffraction [23 that 
LFO orders above room temperature (high for a frus- 
trated magnet) in an AFM state with M;^ = 0, and at 
a lower temperature undergoes a reorientation transition 
to the A' phase inducing a 7^ 0. Our calculations do 
support such a picture, however, we emphasize that this 
wFM state (A' can be thought of as a linear combination 
of Ai and A2) is far more unusual than what could have 
been anticipated. 

The electronic ground state. In the PE phase the 
crystal field at the TM site has a D^h trigonal point sym- 
metry, which splits atomic Zd levels into three sets of 
states as shown in the insets of Figj2]3. In the order of in- 
creasing energy these are: e' doublet composed of dxz and 
dyz states, e" doublet composed of d^y and (ia,2_ ^2 states, 
and a'l singlet that corresponds to ^3^2 state. It is conve- 
nient to adopt this notation also in the low-temperature 
phase, even though in this case the crystal field has lower 
symmetry and the doublets split into singlets. The den- 
sity of states (DOS) plots calculated for LFO and LMO 
in the FE phase are shown in Fig. LFO is a charge- 
transfer insulator with the conduction band formed by 
minority Fe 3(i states and the valence band composed 
of O 2p states, below which are the filled majority Fe 
3(i bands. In the case of LMO majority 3(i bands are 
partially filled with occupied e' and e" . Due to smaller 
nuclear charge of Mn these bands lie at higher energies 
as compared to LFO and they overlap with O 2p states 
forming the valence band. At the conduction band min- 
imum we have the empty majority a'^ states while mi- 
nority 3(i levels lie at somehow higher energies. The im- 
portance of these differences will be made clear when 
discussing the magnetic interactions. 

Symmetric exchange. There are two important 
symmetric exchange interactions. The first is a strong. 
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AFM superexchange interaction between in-plane nn 
spins, Jnn- Its magnitude is much larger for LFO com- 
pared with LMO, suggesting a substantially larger mag- 
netic ordering temperature for ferrites (within mean field 
theory the calculated Currie- Weiss temperature for LFO 
in the FE phase is Qcw = 1525 K, while in LMO 
6cH- = 274 K [34 ). 

The second is a weak, super-superexchange interaction, 
Jc, which couples consecutive spin planes via a TM-0- 
Lu-O-TM exchange pathway. In the PE structure, a spin 
in one layer is connected to three spins in a consecutive 
layer. Each of these degenerate spin- spin interactions has 
two equivalent exchange pathways. We find that the in- 
terlayer exchange is AFM for LFO but FM for LMO. Al- 
though the strength of this interaction is relatively weak, 
this sign difference turns out to be key. 

In the PE structure symmetry implies that the relative 
orientation of the spins in consecutive spin planes is arbi- 
trary. The trimer distortion, however, splits the three de- 
generate interactions into: a single J^^ interaction, medi- 
ated by two equivalent TM-O-Lul-O-TM exchange path- 
ways, and two J^^ interactions, where each interaction is 
mediated by a TM-O-Lul-O-TM and a TM-0-Lu2-0- 
TM exchange pathway. This remarkably introduces an 
extra contribution to the energy 

^f,fra = 2AJ,COs(xi-X2) (3) 

where AJc = J^^ — J^^, for which the sign is key in 
determining the magnetic configuration type: A- type 
(X2=Xi+7r) for AJc > or B-type (x2=Xi) for AJ^ < 0. 

A simple structural analysis shows that the super- 
super exchange mediated through the Lu2 ion is always 
weaker than that mediated through the Lui ion, and in- 
deed our calculations show that the magnitude of J^^ 
is always larger than J^^. We therefore see that the 
choice between A-type and B-type in the FE structure 
is in fact determined by the sign of Jc in the PE struc- 
ture. This is important. In ferrites, the AFM nature of 
the interlayer exchange is uniquely determined by the or- 
bital occupancy, it is always AFM and therefore ferrites 
will always prefer A-type magnetic configurations and the 
wFM ground state. (Although the interlayer exchange in 
LMO is FM, which explains why it prefers B-type mag- 
netic configurations, it is not universally so; a discussion 
is given in the Supplement). 

Dzyaloshinskii-Moriya and Anisotropy. Sym- 
metric exchange determines the magnetic configuration 
type. It is the DM interactions and SI A that drives spin 
canting. (Note, our calculations reveal that a dominant 
DM interaction, D^j, exists only between nn spins within 
the triangular planes.) Let us first consider the PE struc- 
ture. Here the DM vector has only a z component, which 
favors spin ordering within the xy plane while the SIA 
tensor is diagonal with only one independent parameter. 
In the FE structure, however, the trimer distortion 
induces both a transverse component of the DM vector. 



D^j^, parallel to the xy plane and an off-diagonal compo- 
nent to the SIA tensor, ^35j . 

These trimer-induced interactions lead to a further en- 
ergy lowering if the spins in each layer, a, cant out of their 
respective xy plane. The canting angle in each layer is 
determined by minimizing Eq. |2] and is found to be 

M-octan(2^J^ ('^-^-^-f-^^°^(^"\ (4) 

oJnn 

for the ^ = domain, which vanishes in the PE phase 
(as Dj^y (X <I>^2 (X Qks 0, where D^^y |D^J|). 

Note that this result is easily generalizable to any ^- 
domain. Next we consider a simple phenomenological 
model that provides a more intuitive interpretation and 
rigorous generalization of these spin model results and 
will clearly show that the trimer-induced exchange inter- 
actions are responsible for the bulk P-M coupling. 

The phenomenological model. The crystallo- 
graphic structure suggests that insight may be gained by 
considering order parameters, OPs, local to each of the 
two bipyramid layers. Within each layer, a, a local mag- 
netic structure can be defined as a combination of 120° 
AFM order in the xy plane, = L^yC'^^'^ ^ and a net 
magnetization along M^. A specific choice of angles 
(XI7X2), which are defined by specific TM ions from ad- 
jacent triangular planes connected by the 2c axis that re- 
mains in the FE phase, describes any possible AFM spin 
pattern shown in Fig. [2^. Additionally there are two ways 
to combine the layer magnetizations: = ^Mf and 
= Mf - Mi. 

Following a similar line of reasoning, we describe the 
total trimer distortion, 0%^^ as the sum of two local ones, 
Qa = Qks^^^'^i where we generalize the suggestion by 
Artyukhin et al., [16 and define as the angle that 
describes the in-plane displacement of the apical oxygen 
that lies directly above the reference spin, Lq,, as shown 
in Fig. [3J1. Note that ^1 = tt - ^2 = ^ is fixed by the 
symmetry of the trimer distortion, therefore specific ^^'s 
are associated with the different trimer domains, Fig.[l^. 

Let us ask the question, when is a (local) trilinear cou- 
pling between Qa, i^a, and either Mz or allowed by 
the symmetry of the AFM-PE reference structure. A 
derivation is provided in the Supplemental but note, only 
invariants that contain (3? = real part), where 

= QiL\ + r]Q2Ll and r] = ± (e.g., pairs like Qa, L*^ 
are required to preserve translational symmetry), can be 
formed. We find two allowed couplings, 

F,^ ^ 5?[X+]M, = (Qi • Li + Q2 • L2)M, (5) 
F^^f ^ 3?[X_]L, = (Qi • Li - Q2 • L2)L, (6) 

with Lq, = Lxy{x cos Xa + ^sinxa) and = 
Q Ksix cos a + ^sin^Q,). The presence of either invari- 
ant indicates that as the system enters the FE phase, the 
trimer distortion induces a canting of the AFM ordered 
spins in each plane, leading to a trimer-induced wFM or 
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wAFM, respectively [36]. Notice that Eq.[5]not only con- 
tains the spin model result, Eq. |4j but also generalizes it 
to any ^-domain. This connection can be made more 
apparent by defining D^^ by Qq, = z x Dq,, then 

[Di-(LiXz)+D2-(L2Xz)]M, (7) 
[Di • (Li X z) - D2 • (L2 X z)]L, (8) 

leading to an induced magnetization, oc (D^ x Lq,). 
We interpret D^^ as an effective transverse DM vector, 
D^^, for each layer that is induced by the local trimer 
distortion (this can be derived, see Supplemental). 

We see that writing in terms of local OPs has sev- 
eral conceptual advantages. First, it is immediately clear 
that in the FE phase canting occurs only when there is 
a nonzero projection of an xy-sp'm along the direction of 
the local trimer distortion, i.e., cx Qq, • Lq,. This is 
why there is no canting for the Ai and B2 spin configu- 
rations, as the local spins are perpendicular to the local 
trimer distortions. In the A2 and Bi configurations, how- 
ever, the local xy-sp'ms and trimer distortions are parallel 
(or antiparallel) , where the only difference being that in 
the A2 state this projection has the same sign in adjacent 
layers while in the Bi state the projection changes sign, 
resulting in a nonzero and F^^f , respectively. Since 
an intermediate spin state can be thought of as a linear 
combination of principle spin structures, any one with a 
nonzero A2 (Bi) component will exhibit wFM (wAFM). 

Implications of the trilinear coupling. The trilin- 
ear coupling of Eq. |5]is quite remarkable. It implies that 
in ferrites the trimer distortion not only induces a polar- 
ization, Eq. [1] but also mediates a non-trivial bulk P-M 
coupling. To make this clear, we consider a thought ex- 
periment in which an applied electric field, can switch 
P to any one of the three trimer domains with — P (more 
in the Discussion). Let the system be initially in the 
domain. Fig. |3]3, with P = +P and Qa • Lq, = +1, apply 
E and consider two cases: 

{a+ a~}. In a proper FE like PbTiOa the struc- 
ture of the +P domain is related to the — P domain by 
a reversal in the direction of the polar distortions w.r.t 
the PE structure, e.g, the Ti^+ ion moving from up to 
down. The analogous situation in the hexa systems corre- 
sponds to a structural change from a 2-up/l-down buck- 
ling and tilting 'out' of the R-planes and bypyramids, 
respectively, to a l-up/2-down and tilting 'in', while re- 
maining in the same distinct domain, e.g., a. This cor- 
responds to switching P via rotating ^ by tt. In this a~ 
domain Qa • L^ = —1, Fig. [Sj), and therefore, because of 
Eq. |5j either has to rotate 180° or the small canting 
angle has to change sign. It is not unreasonable to expect 
the latter to be more favorable, leading to a reversal of 
M^. This is the scenario previously discussed in Fig. [l] 

{a+ /3~}. The improper nature of ferroelectric- 
ity, however, offers an even more interesting possibility 
in that there exists three distinct and accessible domains 
(a, and 7). As an example let P switch via rotating ^ 



P domain ^=7r/3 




Initial configuration 
b a+ domain ^=0 



t 



Immediately after switching Q to 
C CX" domain ^=71; 



P=+P (P= cos [3^]) 
Q= +Q (P ~ Q') 
U= +L (Q.L~ +1 ) 
M= +M (M ~Q.L ) 




P domain ^=7t/3 





^Q.L~-l 

P=-P(P^ cos [3^]) 
Q= -Q (P ~ Q') 
L=+L(Q.L~-1) 
-M (M ~Q.L ) 



A long time after switching to p domain ^=7r/3 
Two possible equilibrium magnetic domains 
C L switches by 60° 



^0<Q.L <+l 

.". out of equilibrium 
P=+P (P= cos [3^]) 

Q= +Q (P ~ Q') 




L switches by 120° 




^Q.L~+1 

P=-P(P^ cos [3^]) 
Q= -Q (P ~ Q') 
L=-L(Q.L~+1 ) 
M= +M (M ~Q.L ) 
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P=-P (P^cos [3^]) 
Q=-Q (P~Q') 
L=+L(Q.L~-1) 
-M (M ~Q.L ) 



FIG. 3. Bulk cross-coupling of polarization and mag- 
netization in hexagonal ferrites. Thought experiment to 
elucidate the cross-coupling of ferroelectricity and magnetism 
mediated by the trimer distortion for the predicted A2 mag- 
netic ground state, a, Definition of local trimer angles ^1 
and ^2, b, Initial equilibrium domain = {+P,+Q, H-L^y, 
+M;2}. Immediately after switching c, to the a~ domain, 



Lxy still in ground state, therefore -\-Mz 



-Mz\ a domain 



= {-P,-Q, +L^^, -M4, d, P domain, y is not in a ground 
state and must rotate either, e, 60°, therefore -\-Mz +M2; 
/3" domain = {-P,-Q, -L^y, +M4, or f, 120°, therefore 
^ -M^; /3- domain = {-PrQ, +Lxy, -M^. 



by 7r/3 and consider the configuration immediately after, 
Fig. [Sp. In this (3~ domain < Qa • Lq, < +1, implying 
the system is not in equilibrium, and therefore La must 
rotate by either 7r/3, Fig. [3^, or -27r/3, Fig. [sf. In the 
former case Qa • = +1 as in the initial a~ configu- 
ration, therefore is not reversed, while in the latter 
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Qa • = — 1 and switches 180°. 

Discussion. In this paper we have discussed an in- 
triguing family of improper ferroelectrics where we have 
shown for the first time that a non-polar trimer structural 
distortion not only induces an electrical polarization but 
also induces weak-ferromagnetism. It is a universal fea- 
ture of these hexa ferrites systems that they stabilize in 
one of the A- type states, all of which (except for one, the 
Ai state) allow for a weak ferromagnetic moment along 
z. Due to a competition between spin canting and in- 
plane anisotropy, the magnetic ground state is the A2 
magnetic phase. It is important to note that exchange 
interactions in these ferrites are much larger compared 
with the manganites, which indicates the possibility for 
a much higher magnetic transition temperature and an 
above room temperature multiferroic, the observation of 
which has been recently suggested by [23 . 

Here, however, we discussed a result even more ex- 
citing - the trimer distortion mediating an intrinsic bulk 
trilinear-coupling of the polarization, magnetization, and 
antiferromagnetic order. Due to the improper nature 
of ferroelectricity, the topologically protected trimer do- 
mains have a particular phase relationship that leads to 
an appearance of structural vortices, which in turn induce 
magnetic vortices [16]. In these ferrites we additionally 
have the nontrivial bulk magnetoelectric effect, and it is 
natural to ask what type of domains would one expect 
to see here? Here we dsicuss two likely domain patterns. 
Fig. [4j but by no means is this meant to be exhaustive. 
Both assume that the only AFM domain walls are those 
clamped to the trimer domain walls [16] . 

The first domain configuration. Fig. |4^, is such that 
at trimer domain walls differing by = 7r/3, e.g., 
at the a+ and /3~ domain wall, the AFM spins rotate 
^Xa = k/3|. In this case the direction of the magne- 
tization remains the same across the domain wall as we 
previously explained in Fig. [3^. Note that this seem- 
ingly lower energy AFM switching pathway results in a 
homogeneous magnetization across the entire material, 
which has to be unstable towards the formation of ferro- 
magnetic domains. But what kind of domains? In these 
hexa systems iM^ domains result from the reorientation 
of the Lixy vector by an angle tt. They can occur within 
the bulk of a trimer domain, i.e., a free AFM domain, 
but there is an energy cost to form the domain wall. 

An alternative path to minimize the total energy of 
the system is considered in Fig. [4]3, where at trimer do- 
main walls that differ by = 7r/3 the AFM spins now 
rotate Axa = |27r/3|. In this case, the magnetization 
direction reverses with the polarization as discussed in 
Fig. |3]F. Therefore, even though this domain configura- 
tion at first appears less likely than Fig [4^, it provides 
an avenue for the system to minimize the magnetostatic 
energy without having to introduce free domains. 

Finally, the hexa ferrites provide an intriguing route 
to a bulk magnetoelectric effect and electric-field control 




FIG. 4. Clamping of trimer and magnetic domains and 
electric field response. At the lowest energy trimer domain 
wall ^ changes by 7r/3, which will rotate the AFM order pa- 
rameter Li and L2 by either a, 7r/3 and —tt/S or b, — 27r/3 
and 271/ 3. In the former case remains constant while in 
the latter reverese at each trimer domain wall, i.e., fer- 
roelectric, ferromagnetic, and trimer domains are mutually 
clamped, c, The schematic diagram displaying the electric- 
field switching of ferromagnetic domains (note, the figure is 
inspired by Ref . pT| ) 

over the direction of the magnetization. In Fig. ^ 
we sketch the expected response of the domains to 
electric-field poling. (Note that as the domain walls 
are topologically protected, however, -P domains do 
not disappear completely [TT| [16^.) In this process the 
positive electric field, E, e.g., choses the (+P,+M) state 
and therefore reversing of the direction of electric field 
will not only switch the direction of polarization, but 
also reverses the direction of magnetization. 
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